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ABSTRACT 

We propose a new heating mechanism in magnetar crusts. Magnetars' crustal magnetic fields are much 
stronger than their surface fields; therefore, magnetic pressure partially supports the crust against gravity. The 
crust loses magnetic pressure support as the field decays and must compensate by increasing the electron degen- 
eracy pressure; the accompanying increase in the electron Fermi energy induces nonequilibrium, exothermic 
electron captures. The total heat released via field-decay electron captures is comparable to the total mag- 
netic energy in the crust. Thus, field-decay electron captures are an important, if not the primary, mechanism 
powering magnetars' soft X-ray emission. 

Subject headings: nuclear reactions, magnetic fields — stars: neutron 



L EVITRODUCTION 

Soft gamma repeaters (SGRs) and anomalous X-ray pul- 
sars (AXPs; for reviews, see IWoods & Thomps"onl 120061: 
iKaspil I2007t iMereghettil 120081) are the two known mani- 
festations of magnetars: neutron stars with ultrastrong ex- 
terior magnetic fields B > lO'"* G and eve n stronger inte- 
rior, primarily toroidal fields (Duncan & Thompson 119921: 
iPaczvriskil 19921; [Thompson & Duncan 1993) . As of this writ- 
ing, all magnetars have characteristic ages Tchar ^10^ yr, and 
their persistent soft X-ray lumin osities L ^ 10^^ ergs"' (al- 
though some are transients, e.g.. llbrahim et al. 2004), which 
are likely of thermal origin, are systematically higher by 
factors of > 100 than the power ava ilable from rota tion 
jvan Paradiis et al.lll995h. cooling (e.g.. iPage et al.ll2006l) . or 
accretion (i Hulleman etani2000l: iKaplan et al.ll200lir Taken 
together, these observations suggest that magnetars' magnetic 
fields de cay and thereby induc e heating in the stellar interi- 
ors (e.g .,lT hompson & Dunca n' 1996; Hevl & Kulk a rni,,1 998; 
IColpi e t al. 2000; Pons & Geppert 2007; Pons et al.l l2007l)" 

Sustaining the persistent thermal emission requires an inte- 
rior energy reservoir 



£inf > 3 X lO"* 



"^char 



105 yry V 1035 ergs- 



erg. 



(1) 



Comparing this to the magnetar' s total internal magnetic en- 
ergy Eb ^ (B^/8Tr){47rR^/3) gives B > 1 x lO'^ G. However, 
these values are likely flagrant underestimations: First, a mag- 
netar emits most of its therm al energy as neutrinos. Indeed, 
IKaminker etaH (l2006l l2009l) argue that, to avoid severe neu- 
trino losses, the heat source powering the thermal emission 
must be located at or near the outer crust, i.e., within the mag- 
netar's outermost « 100 m. If the thermal emission originates 
from magnetic energy contained only within the outer crust, 
then B ^ 10'* G. Second, a magnetar's magnetic field powers 
other phenomena as well, such as giant flares; the 2004 De- 
cember 27 giant flare from SGR 1 806-20 released 10^* erg, 
which implies B ~ 10'* G (assuming giant flares are recurre nt; 
IHurlev et allllOOl ISteUa et alJl2005tlTerasawa"etani2005h . 

That the magnetic field's strength is much greater in the 
crust than at the surface implies that a substantial B gradi- 



ent exists. Therefore, if the crustal field geometry is primar- 
ily nonradial, magnetic pressure partially supports the crust 
against gravity. We argue in ^that the crust loses magnetic 
pressure support as B decays and compensates with an in- 
crease in electron degeneracy pressure. The accompanying in- 
crease in the electron Fermi energy induces exothermic elec- 
tron captures that heat the crust. We calculate the total energy 
released via field-decay electron captures in |3]and compare 
it to Joule heating in ^ We discuss our results and speculate 
on other consequences of field-decay electron captures in ^ 



2. FIELD-DECAY ELECTRON CAPTURES 

We construct a simple model of the hydrostatic struc- 
ture of a neutron star's strongly magnetized outer crust to 
illustrate the process by which magnetic field decay in- 
duces electron captures. We ignore the inner crust both 
for simplicity and because most of the thermal energy re- 
leased there escapes as neutrinos. Consider a neutron 
star of mass M and radius R. The outer cru st's scale 
height (rs 100 m, e.g., IShapiro & Teukolskvlll983h is much 
less than R; therefore, we adopt a plane-parallel geom- 
etry and set the gravitational acceleration g = GM/R^ to 
be constant throughout the layer. Furthermore, we make 
three assumptions regarding the crustal magnetic field: (1) 
Its nonradial (e.g., toroidal) component is large relative to 
its radial component, (2) the field strength in the crust 
is greater than that at the stellar surface, and (3) the 
field decays with time. These assumptions are likely 
satisfied in magnetars ([ Thompson & Duncan 1993, 1995, 
'1996, 2001; Braithwaite & Spruit 2006; Geppert et al. 200a 
Pere z-Azorm et al. 2006a b; Braithwaite 2009; Reiseneggeg 
120091) . 

The equation of hydrostatic equilibrium in the electron 
magnetohydrodynamic limit is 



Pg-V^-V(^)+^(B.V)B = 0, 



(2) 



Hubble Fellow 



where p is the density, P is the pressure, and B is the mag- 
netic field. We ignore the magnetic field's affect on P (see 
Appendix). Setting the radial component of B equal to zero 
by assumption (1) for simplicity, the equation of vertical hy- 



2 



drostatic equilibrium becomes 



Pg-- 



— (p+ — 



(3) 



Invoking the plane-parallel geometry and rewriting Equation 
(O in terms of the Lagrangian coordinate S, the column depth 
as measured from the surface (such that t/S = -pdr). 



(4) 



Finally, integrating Equation (|4|i from the surface, we get 



^E = P(I],f)+— [B'(S,f)-B'(0,f)], 
ovr 



(5) 



where we have set f = at the surface, where E = 0, and 
included the explicit time and spatial dependences of P and 
B. Equation (|5]l shows that electron degeneracy pressure and 
the magnetic field together support the column of matter E 
against gravity. 

Consider a matter element at column depth E containing 
a nucleus of mass M(A,Z), where A is the mass number and 
Z is the proton number, as B decays with time. The column 
gradually loses magnetic pressure support; to compensate, P, 
and thereby the electron Fermi energy E-p, increases with time. 
When M(A,Z)c^ + Ep > M{A,Z- l)c^, electi-on capture onto 
the nucleus becomes energetically favorable. 

The nuclear species at a given E is a function of Ep (or 
equivalently, P; e.g., Shapiro & Teukolskv 1983), so even 
when the crust is initially at equilibrium, electron captures 
will occur: Any nucleus that is stable at a given Ep is unsta- 
ble to capture at a sufficiently larger Ep. The loss of magnetic 
pressure support due to field decay forces nuclei out of their 
initial equilibrium and into a new equilibrium configuration. 

Electron captures can be exothermic in two mutually inclu- 
sive ways (e.g., Gupta et al. 2007): (1) An electron captures 
into an exited state of the daughter nucleus if, for example, the 
daughter nucleus's ground state is forbidden. The daughter 
nucleus then radiatively deexcites and thereby heats the crust. 
(2) If the parent nucleus is even-even, then M{A,Z- 1) > 
M(A,Z-2) due to the nuclear pairing energy, and a second, 
nonequilibrium electron capture immediately ensues and re- 
leases heat. 

This process resembles the well-known deep crustal heating 
in accreting neutron stars (e.g. .'Bro wn et al.ll998tlGupta et alj 
12007, 2008; Haensel & Zdunik 2008]): In that context, P alone 
supports the matter column against gravity, so ^E = P. The 
accretion-driven advection of a matter element to higher pres- 
sures, and thereby higher Ep, induces electron captures that 
heat the crust. In magnetars, however, it is not accretion, but 
magnetic field decay, that induces electron captures^. 

3. CRUSTAL HEATING VIA FIELD-DECAY ELECTRON CAPTURES 

We now estimate the total thermal energy released in the 
outer crust via field-decay electron captures. Consider a neu- 
tron star crust with an initial magnetic field B, and that con- 
tains an ion species with an electron capture threshold energy 
^F.ec, such that electrons capture when Ep > /sp.ec- After a 

- The effective gravitational acceleration increase during spin-down of a 
rapidly rotating neutron star raises Ej: of a matter element and thus induces 
electron captures as well, although it is unimportant for SGRs and AXPs 
because they rotate slowly. This effect is similar to rotochemical heating in 
neutron star cores tReisenegger..l995,) . 
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Figure 1. Schematic diagram of the electron Fermi energy Ep as a function 
of column depth E during crustal magnetic tield decay. Electron degeneracy 
pressure P and magnetic pressure B^/Sir together support the column against 
gravity. The column loses magnetic pressure support as B decays from Bi 
to Bf, and P and Ef concomitantly increase. When > Ep ^c, the electron 
capture threshold, exothermic captures ensue and heat the crust. The matter 
column within which field-decay electron captures occur is Sec,i — Sec.f- 

period of time, the crustal magnetic field decays (via unspeci- 
fied means) to Z?f <C Bi. For relativistic, degenerate electrons, 
P cx Ep; we define Pec to be the pressure at which Ep = /ip.ec- 
We now determine the column of matter within which elec- 
trons capture as the field decays from Bi to Bf. The vertical 
hydrostatic equilibrium equations for columns of matter at the 
electron capture pressure Pgc before and after field decay are 
(Equation ®) 



gEeci =P.c+^ [B?(Sec,i) - Bfm , 

^Eecf = P.c+^ [B?(Sec,f) - Bjm , 
OTT 



(6) 
(7) 



respectively, where Eec.i and Eec,f are the electron capture col- 
umn depths. Subtracting Equation (|7]) from (|6]l and noting that 
Bi(Sec,i) > Bi(0), Bf(Eec,f), and Bf(0) by assumptions (2) and 
(3) gives the total column of matter available as nuclear fuel. 
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(8) 



See Figure[T]for a schematic diagram. 

Assuming spherical symmetry, the total energy released via 
electron captures 



; 47r/;2(Ee 



mu 



(9) 



where Q is the energy per nucleon released via electron cap- 
ture and niu is the mean nucleon mass. In general, a nu- 
cleus will undergo several electron captures in the outer crust 
as the field decays, but for most rea ctions 2 < 0.1 MeV 
dGupta et all 120071: iHaensel & Zdunikll2008i) . Captures that 
release the most heat occur near neutron drip (where P w 
1 X 10^ ° dynecm"^), i.e., those that induce neutron emissions; 
iGupta et al., (,2008.) find Q ~ 1 MeV. From Equations ^ and 
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(|9|l, and normalizing to conditions at neutron drip, 

2 



; 3 X 10^ 



1016 G 



1 MeV / 



R 



10 km 



M 



1.4M0 
(10) 



4. COMPARISON TO JOULE HEATING 

The usual explanation for magnetars' persistent thermal 
emiss io n is Joule hea ti ng in th e crust (e.g.. IMi ralles et al. 
T99l lAguilera et al.1 EOOl lUrpin & Koneiikovi i2008i: 
Pons et al .I l2009h : The crust's electrical resistivity is nonzero, 
so currents generating the internal magnetic field decay in 
time, converting some of the magnetic energy to thermal 
energy and thereby heating the crust. In this section, we 
compare heating via field-decay electron captures to Joule 
heating. 

From Equation (|8]l and the expression for the capture pres- 
sure scale height for relativistic, degenerate electrons //ec = 
Pec/pg = £^F.ec/4/ie'Wu^, whcrc /ig is the mean molecular 
weight per electron. Equation (|9]l becomes 



i 4ttR^H„ 



4/iee 



■^F.ec 



(11) 



The term 47r/?^//ec(Bf /Stt) w Eb, the total magnetic energy 
contained in a crust of scale height i/ec; the ratio of the total 
nuclear energy released via field-decay electron captures to 
the total magnetic energy 
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l^ = 0.5f^ 

£^F.ec V3.2 



Q \ ( £F.ec y 

1 MeV ) V27MeV / 



(12) 



Nuclear heating and Joule heating are comparable if all of the 
magnetic energy converts to thermal energy. However, un- 
like electron captures, the crustal magnetic field powers phe- 
nomena other than the thermal emission, such as the hard (~ 
100 keV) nonthernial X-ray emission (e.g..lKuiper et al.l2004t 

■ m 



iMolkov et al.ll2004tlRevnivtsev et al.l20b4ir X-ray bursts, and 
giant flares. The fraction of magnetic energy that contributes 
to the thermal emission is unknown. Therefore, we conclude 
that heating via field-decay electron captures is a consider- 
able, if not the primary, source of the persistent soft X-ray 
luminosity from magnetars. 

However, we stress that equations ( [TOlfTSb are upper limits 
with respect to our assumed magnetic field configuration. If 
the field were instead predominantly radial or force-free, for 
example, field decay would induce httle nuclear heating. 

5. DISCUSSION AND FUTURE WORK 

In this work, we identified a new, strong heat source in mag- 
netar crusts. Magnetic pressure partially supports the crust 
against gravity. The electron degeneracy pressure increases 
as the crust loses magnetic pressure support during field de- 
cay, and the concomitant increase in induces exothermic 
electron captures. The total heat released via field-decay elec- 
tron captures is of order the crust's total magnetic energy. 
Field-decay electron captures are important to neutron stars' 
magneto-thermal evolution and should be included in future 
calculations. At a minimum, field-decay electron captures 
supplement Joule heating and lessen the interior magnetic 
field strength required to power magnetar phenomena; how- 
ever, since only a fraction of the crust's magnetic energy, but 
essentially all of the nuclear energy, converts to thermal en- 
ergy, field-decay electron captures may in fact dominate mag- 
netars' persistent soft X-ray emission. 



We now outline other potential ramifications of this work: 
_[ (1) Although we neglected both the inner crust and outer- 
most envelope for simplicity, field-decay-induced nuclear re- 
actions can occur there as well. In addition to electron cap- 
tures in the inner crust, the loss of magnetic pressure support 
may trigger light-element fusion reactions in the envelope or 
pycnonuclear reactions in the inner crust, since p increases at 
a fixed S during field decay. 

(2) The crust's elec trical resistivity decreases with p and in- 
creases with T (e.g., Yakovlev & Urpin'lQSO"), so magnetic 
field dec ay via Hall-drift-enhanced Ohmic dissipation (e.g., 
iGoldreicii & Reiseneg ger 1992) is most efficient at low densi- 
ties and high temperatures. Magnetic pressure support lowers 
p at a given S relative to the nonmagnetic case, and field- 
decay electron captures heat the crust; both effects accelerate 
Ohmic dissipation. 

(3) Although it is the dominant mechanism, field decay is 
not required to induce electron captures. For example, rear- 
rangement of the crustal magnetic field via Hall drift, which is 
nondissipative, can induce captures: Altering magnetic pres- 
sure gradients alone can either raise or lower the electron pres- 
sure, and thereby E-p, of a matter parcel past an electron cap- 
ture threshold. Raising /ip past a capture threshold triggers 
exothermic electron captures, but lowering Ep past a thresh- 
old does nothing, since exothermic captures are effectively 
irreversible. This may further enhance the heating rate, and 
hence the Ohmic dissipation rate, in the crust. 

(4) Asymmetries in the internal magnetic field relative to 
the stellar rotation axis produce distortions in a young, rapidly 
rotating magnetar that may generate gravitational waves. A 
nascent misalignment between the rotation and magnetic axes 
is secularly unstable a nd grows with time (e.g.. lCutTe3l2002l ; 
iDairOsso et al.l 120091) . If the distorting magnetic field is 
predominantly nonradial and within the crust, the resulting 
field evolution causes electron captures, which can gener- 
ate asymmetric densi ty variations that intensify gravitational 
wave emission (see Bildstenlll998t lUshomirskv etaLl llOOOl 
who proposed a similar mechanism for weakly magnetic, ac- 
creting neutron stars). 

(5) We neglected the crust's elasticity in the hydrostatic 
equilibrium equation (|2|i. Elasticity partially compensates lo- 
cal magnetic stresses and thereby delays electron captures; 
conversely, sudden pressure increases that offset fracture- 
induced elasticity losses during catastrophic events like SGR 
bursts and giant flares may trigger electron captures. While 
unlikely to contribute to the burst's emission itself, captures 
triggered by crust frac ture may contribute to the subsequent 
transient cooling (e.g.. lLyubarsky et al.ll2002l) . 
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Figure 2. Electron Fermi energy as a function of pressure P = Pe + , 
where Pc is tlie electron degeneracy pressure and Pc is the Coulomb correc- 
tion, for different magnetic field strengths B. A magnetic field decreases at 
fixed P relative to the nonmagnetic case; Ef increases as B decays, inducing 
electron captures. The kink in the B > plots indicates the Ef valu e below 
which electrons occupy only the ground Landau level (Equation |A2)- 

APPENDIX 

ELECTRON CAPTURES VIA LANDAU QUANTIZATION LOSS 

A magnetic field quantizes electron motion perpendicular to 
B into Landau orbitals and thereby lowers the electron degen- 
eracy (e.g., Harding & Lai 2006). The electron degeneracy 
pressure oc n^E-p, where is the electron number density; 
for a given P^, B increases «e and lowers E-p relative to the non- 
magnetic case. Therefore, Ep increases as the magnetic field 
decays, inducing electron captures. The Coulomb correction 

4/3 

to the pressure Pc oc -nj " negligibly mitigates this effect. 

Landau quantization affects Ep only when P <C B^/Stt (Fig- 
ure|2]i. The electron energy eigenstates 



J J 7 4 



\+2v- 



1/2 



(Al) 



where p±^ is the electron momentum perpendicular to B, v 
is the quantum number, Bcrit = m\c^ /Tie = 4.414 x 10'^ G is 
the quantum critical field strength, and me, c, h, and e are 
the usual fundamental constants. The magnetic field has the 
largest effect on P when all electrons are in the ground Landau 
level v = Q. From Equation dAlb . this occurs when 



Ep<m,c' ( 1+2 — 



1/2 



10.9 



B 



101'' G 



1/2 



MeV. (A2) 



From Equation iA2\ and the relation Pg ~ £'p/127r^(7zc)^ 
for relativistic, degenerate electrons at Z? = (e.g., 
[Shapi ro & Teuko lskv 1983), Landau quantization affects P 
when 

^<|^ = 6xl0-^ (A3) 

where a = /lie is the fine structure constant. 

Unlike captures via magnetic pressure support loss (©, 
captures via Landau quantization loss always occur during 



magnetic field decay, i.e., they occur independent of the field's 
geometry. However, Equation (IA3b shows that Landau quan- 
tization's effect on electron captures is small relative to that 
of magnetic pressure support if B is primarily nonradial. 
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